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We show that, due to the weak ferromagnetism of I^-zSrajCuCU, an external magnetic field 
leads to a dimensional crossover 2D — > 3D for the in- plane transport. The crossover results in an 
increase of the hole's localization length and hence in a dramatic negative magnetoresistance in the 
variable range hopping regime. This mechanism quantitatively explains puzzling experimental data 
on the negative magnetoresistance in the Neel phase of I^-xSrajCuO.!. 



I. INTRODUCTION 

The physics of the high-temperature superconduct- 
ing oxides is determined by the interplay between the 
charge and spin degrees of freedom, ultimately respon- 
sible for the superconductivity itself. A variety of inter- 
esting phenomena exists already at low doping when the 
oxide layers are insulating. In La2- K Sr 2 .Cu04 (LSCO), 
the insulating (spin-glass) region corresponds to doping 
x < 0.055, with incommensurate magnetism which ex- 
ists down to the boundary with the antifcrromagnctic 
phase (at x = 0.02), and even inside the Neel region 
(x < 0.02) pi A popular point of view favors an expla- 
nation of the incommensurate magnetism based on the 
tendency of the holes to form "stripes."— However, experi- 
mental data on variable range hopping (VRH) (see the re- 
view Ref. [H) , unambiguously indicate localization of holes 
for x < 0.055 and therefore support an approach based on 
a purely magnetic scenario, where a spiral distortion of 
the spin background is generated by localized holes. The 
corresponding theory explains quantitatively the variety 
of magnetic and transport data in LSC04£&£&£ 

Magnetic phenomena in the low-doping region re- 
flect, in addition to the Heiscnbcrg exchange, the pres- 
ence of anisotropics in the spin-spin interactions, such 
as Dzyaloshinsky-Moriya (DM) and XY terms. In the 
present paper we consider the Neel phase, x < 0.02. 
In this phase the anisotropics confine the spins to the 
(ab) plane and fix the direction of the Neel vector to 
the 6-orthorhombic axis. Moreover, the DM interac- 
tion induces a small out-of-plane spin component that 
is ferromagnetic in the plane (weak ferromagnetism) but 
staggered in the out-of-plane c-dircction. This com- 
ponent can be easily influenced by an external mag- 
netic field applied in different directions, as it has been 
recently addressed both experimentall y 10 ' 11 ' 12 ' 13 ' 14 and 
theoreticallyi 15 ' 16 For example, a perpendicular field (H j 
c) can cause an alignment of the out-of-plane moments 
via a spin- flop transition at a critical field Hf , determined 
by the competition between the DM and inter-layer 
Heisenberg exchange (typically Hf ps 5 — 7 T) i n ' 13 ' 14 



Perhaps most intriguingly, the in-plane resistivity (along 
with the c-axis resistivity) decreases by as much as 50% 
across such a transition. 10 ' 11 The magnitude of the mag- 
netoresistance (MR) shows a rapid increase only below 
ps 50 Kii where LSCO exhibits VRH conduction^ 
This implies that the MR is accumulated mostly in tran- 
sitions between localized states. Therefore it is very nat- 
ural to assume that the large negative MR is due to an 
increase of the hole's localization length as it was sug- 
gested in the first experimental paper ^ From theoretical 
viewpoint the problem is why the localization length in- 
creases at the spin flop transition. The first model for the 
localization length increase, invoking a three-dimensional 
(3D) VRH mechanism, was proposed in Ref. [ID. How- 
ever, it is clear now that except for ultra-low tempera- 
tures (that we estimate to be below ~ 50mK), the VRH 
conduction at zero magnetic field is dominated by two- 
dimensional (2D) physics* 3 -^ Because of this the 3D pic- 
ture is not able to describe the most recent and detailed 
MR data, as we discuss below. Experiments are per- 
formed typically in the temperature range of a few Kelvin 
and higher where the out-of-plane resistivity anisotropy 
is large p c j p a b ~ 10 2 — 10 3 *ii While we ultimately expect 
that at T — > VRH will become 3D, in the temperature 
range of experimental interest the 2D mechanism is the 
relevant one, as is clear from the analysis of the 2D-3D 
crossover temperature and the fits of the hopping con- 
ductivity presented in the next section. 

In the present work we demonstrate that the large MR 
arises from a change of the effective dimensionality of the 
VRH mechanism with applied field. We support our con- 
clusions by detailed comparison with recent experiments 
on magnetotransport which can be described by our the- 
ory with excellent accuracy. The main idea of the present 
work is that a dimensional crossover (2D — > 3D) occurs at 
the spin flop, and this is conceptually and quantitatively 
different from the 3D picture of Ref. [l8|. In particular in 
our approach the increase of the MR (and the localiza- 
tion length) is not simply due to the change of the out-of- 
plane effective mass as in Ref. [ID, but rather arises from 
a change in the shape of the (localized) wave-functions 
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across the spin-flop. In the temperature regime that we 
keep in mind, IK and higher, the change of the out-of- 
plane effective mass is a small, secondary effect (which 
can manifest itself only at ultra-low temperatures where 
the full 3D VRH mechanism is responsible for transport). 
We show that the alignment of the weak ferromagnetic 
moments in neighboring planes with the field allows the 
inter-layer hopping of localized holes, which in turn leads 
to an increase of the hole's in-plane hopping probability 
and thus negative MR. The presence of an inter-layer 
hopping channel across the spin-flop was already identi- 
fied in Ref. [lH however our analysis differs in the effects 
this additional channel can produce in VRH conduction. 
By investigating the evolution of the hole bound state as 
a function of magnetic field and temperature, we find that 
in various regimes different numbers of layers are involved 
in transport. In the experimentally relevant temperature 
range the hopping turns out to be quasi- two-dimensional, 
leading to a negative MR in very good agreement with 
the most recent experiments ^ 1 ! 12 

The paper is organized as follows. In Section II we 
analyze the effect of the magnetic field on the dispersion 
of the localized holes, through the inter-layer hopping. In 
Section III we present a detailed analysis of the change of 
the hole's wave-function, due to the modified dispersion. 
In Sections IV and V we then use the wave-functions 
to calculate the magnetoresistance for out-of-plane and 
in-plane magnetic fields, and compare with experiment. 
Section VI contains our conclusions. 



II. INTER-PLANE HOPPING AND SPIN-FLOP 
TRANSITION FOR MAGNETIC FIELD H \\ c 

First we briefly summarize previous results related to 
the structure of the hole's 2D bound state at zero field. 
As a starting point, we consider the hole dynamics in the 
antifcrromagnetic background within the framework of 
the t—t' — t" — J models In the absence of the Coulomb 
potential V(r) of the Sr ion, a hole resides, in momen- 
tum space, near the nodal points (±7r/2, ±7r/2) and has 
dispersion €k ~ ^-fcj + where /3i f=a /?2 = /3 = 

m^ 1 w 2 J is the inverse 2D effective mass appropriate 

for LSCO& (to|| as 2m e in absolute units). We measure 
energies in units of J = I30meV, and the lattice spacing 
is set to unity. Due to V(r), the hole is localized, and its 
wave function has the form ip{ r ) ~ e _K ° r , corresponding 
to binding energy eo = (3k^/2. Here the inverse (2D) 
localization radius for LSCO is ~ 0.3 — 0.4^ giving 
eo ~ 1 OmcV. On a perfect square lattice the bound state 
is four-fold degenerate: the hole can reside on either up 
or down sub-lattices (pscudospin), and it can reside in 
either of the two pockets {n/2, ±7r/2) (flavor). The or- 
thorhombic distortion lifts the flavor degeneracy due to 
the presence of diagonal next-nearest neighbor hopping 
t 1 . Hence, the holes occupy only the pocket (tt/2, —ir/2)& 
This is also consistent with the fact that the spin struc- 



ture becomes incommensurate along the b orthorhombic 
direction, as seen in neutron scattering for x < 0.055J^ 

Now let us consider the correction to the 2D 
dispersion Se^ arising from the inter-layer hopping 
t±. Without account of correlations we have 5e^ = 
—8t±cos(k x /2)cos(k y /2)cos(k z ), with k-dependent t±, 
t± = t c j(cosk x — cosky) 2 , where t c ~ 50meVJ£ By av- 
eraging over the momentum distribution (k x ,k y ) in the 
2D bound state, we find the effective value of inter-layer 
hopping t± — > £ c Kq/4 ~ ImcV. Since this is a crude, 
order of magnitude estimate, below we will use t± as a 
fitting parameter. 

The t — J model correlations change Se^. First, the 
hopping matrix clement should be replaced by t j_ — ► Zt± , 
where Z as 0.3 is the quasiparticle residue. Second, direct 
hopping is allowed only between spins in the same sub- 
lattice. In LSCO a spin in a given plane (e.g. spin "I" in 
Fig. 1) interacts with four others in the plane above (and 
below) p but at H = (or H < Hf) only out-of-planc 
hopping in the b direction is allowed, because this corre- 
sponds to ferromagnetic ordering of spins in neighboring 
planes (see Fig. I). However, when a magnetic field is ap- 
plied along the c-axis the spins in the next layer reverse 
their signs across the spin-flop transition at Hf, so that 
for H > Hf only hopping in the a direction contributes. 
This is schematically shown in Fig. 1, and reflects in the 
effective dispersion: 

Se u = -AZt ± cos(fc z ) cos (j^A , (i) 
"-" :H <H f , " + " : H > H f , 

where we define the z-direction z \\ c. Since in the occu- 
pied pocket k x as tt/2, k y as —tt/2, Eq. ([I]) reads 

8e" <Hf « 0, H < H f 

5e" >Hf = -iZt±cos(k z ), H > Hf. (2) 

Thus at H < Hf there is no z-dispcrsion, the coherent 
dynamics is purely 2D and the VRH in-plane resistivity 
behaves as p = po exp (To/T) 1 / 3 ^ T is strongly doping 
and sample dependent, and the data at zero field, x = 
0.01 can be fitted by p = 8 x lO^cm and T = 3.6 x 
I0 4 K with astonishing accuracy in the range 4K < T < 
50K, as shown in Fig. [^Inset). 

Observe that a more accurate estimate of the z- 
dynamics below the spin flop can be obtained by expand- 
ing the in-plane dispersion + Se^ around (vr/2, —tt/2) 
and minimizing the resulting quadratic form. This gives 
a non-zero dispersion in the z-direction, but the corre- 
sponding effective mass is huge, Mj_ = f3/{8(Zt±) 2 ] ~ 
I0 4 m e . Consequently we find that at temperatures be- 
low T* - 2- 4 k- 3 (to||/M_l) 3 / 2 T - 50mK, the VRH is 
ultimately three-dimensional, but this regime is irrele- 
vant to present experiments. 
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FIG. 1: (Color online.) Spin structure of LSCO. Red arrows 
correspond to H < Hf, the blue arrows in the middle layer 
show the spin reversal for H > Hf. 



III. EVOLUTION OF THE HOLE BOUND 
STATE ACROSS THE SPIN-FLOP TRANSITION 

To qualitatively understand the effect of the change 
of dispersion @ at H > Hf on the in-plane hole dy- 
namics, we first estimate the change in the hole binding 
energy. From ((2]). after the spin flop the edge of the 
continuum (k z = 0) decreases by 4Ztj_, while the abso- 
lute energy of the bound state, to second order in the 
small parameter Zt±/e.Q <C 1, decreases only by amount 
AE ~ (Zt^) 2 /e <C tj_. The binding energy, which is the 
magnitude of the difference between the absolute energy 
and the continuum limit, changes after the flop as 



e « e - iZt ± + (Zt ± ) 2 /e Q w e - AZt ± . (3) 



Within the VRH picture the conduction is proportional 
to the hole's hopping probability, which decays exponen- 
tially away from the donor site due to the hole localiza- 
tion. Thus, the decrease of the hole binding energy (J3]) 
across the spin flop signals an increase of the localization 
length, and in turn an increase of the VRH conductivity. 
This is our central idea that explains the negative MR. 

To make this argument more quantitative, we need 
to compute the change in the hole's hopping probability 
across the spin-flop transition. Let us enumerate planes 
by the index n and assume for simplicity that Sr produces 
a local potential V(r) that acts only in the plane n = 0. 
For a shallow level the exact form of the local poten- 
tial is not important )2i and for simplicity we will use the 
(^-function approximation: V(r) = —^S(r — ro), where 
r is assumed to be smaller than the localization length, 
ro <C 1/kq. The bound state is described by the wave 
function ipn{r) that depends on both n and r. Before 
the spin flop, H < Hf, it obeys the Schrodingcr equa- 
tion: f-f A r - d n0 fS(r - r )J VnW = E ip n (r), where 



Eq = —£q = — /3kq/2. The solution for r < ro is given by 
Iq(k t), and for r > r it reads 



■>Po = 
H < H 




K r > 1 , 



(4) 



where Io(f) and Ko(r) are the modified Bessel functions 
of the first and second kind respectively. For the inverse 
localization length we obtain kq = 2e r 7 exp (— f}/ (2g)), 
where 7 = 0.577 is Euler's constant. However the exact 
dependence of ko on the parameters of the potential is 
not important, since the energy scale that kq determines, 
the binding energy eo = /3kq/2, was extracted from ex- 
periment in earlier work, eo ~ lOmeV (see Section II). 
For H > Hf the Schrodinger equation becomes 

P \ 

-A r + S n0 V(r) j ?p n - 2Zt ± [lp n+ l + ip n -i] = Eip n , 

(5) 

with V(r) = -H{r - r ), and E = -/3k 2 / 2 - Af *er the 
Fourier transform ipn(j~) = Yip 4>p{ r ) eipn > we obtain 

-|A r ^(r) + y(r)Vo(r) = {E + AZt x cosp)^(r) . (6) 



By solving this equation we find 

/ _ , 1 , r, lrj± I \2 

K/K 

(j> P {r) 



l + 2{Zt ± /e y 

^K (r^/ K 2 -(8Zt ± /(3) cosp) . 



(7) 



Thus, as we have already pointed out before Eq. 
after the spin flop the absolute energy E is shifted only 
in the second order in Zt j_ / eo and hence this shift can be 
neglected. Consequently below we set k — kq. However, 
in contrast to Eq. (U]), the new wave function, 



ipn(r) 



H > H f , 



dp 
2^ 



^ n K ( Ko r^/l - (4Zi ± /e )cosp), 



(8) 



does not have a simple exponential decay. 



IV. MAGNETORESISTANCE ACROSS THE 
SPIN-FLOP TRANSITION FOR FIELD H \\ c 

To evaluate the MR across the spin flop we compute 
now the change of the hole's probability to propagate at 



1 Oh} 1 / 3 
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We 



the (large) VRH distance Rt'- kqRt — ■§ (j?- 
can then identify two large-distance regimes for the wave 
function ([5]): (1) Very large distances, n r ^> eo/(2Zt±), 
and (2) Intermediate large distances, eo/(2ZijJ ^> 
K r > 1. 

In the first regime the integral in (jHJ) can be evalu- 
ated using the saddle-point approximation and the wave 
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FIG. 2: (Color online.) The MR jump \ Ap/p\ (solid, blue line) 
across the spin-flop transition for H\\c. Here Zt±/eo = 0.057, 
T = 3.6 x 10 4 K. The circles are data from Ref. [ll|. Inset: 
Zero field in-plane resistivity fit to 2D VRH form (solid, red 
line) with T = 3.6 x 10 4 K. 



function is spread over many transverse channels, n ~ 
\f AKorZt±/eo 1. The probability density at distance 
r, P(r,tx) = EnhMOI 2 ,k 



P(r, t. 



2r ^/8Tr(ZtJe )n r 
for n r > e Q /(2Zt ± ). 



-2K r- v /l-(4Ztj./eo) 



(9) 



In this regime the wave function has a pure exponen- 
tial decay with a localization length corresponding to the 
shift of the binding energy given in Eq. ([3]). We then find 
that the ratio of conductivities after (H > Hf) and be- 
fore (H < Hf) the flop is 



P(Rr,t±) 



exp 



AZt ± 
3e 



l/3~ 



(10) 

with Vo from Eq.®. This 

- (T„/T)V3 » e /(2Zt x ), 



<7H<H f P(R T ,t± = 0) 

Here P(R T ,t± = 0) = \ip \ 2 
result is valid when kqRt ■ 
which happens in practice when the temperature is below 
IK. In this case UH>H :f I&h<h s ^ 1, i.e. the correspond- 
ing MR is large: \Ap/p = {ph>h s / Ph<h s ) ~ 1| ~ 1, and 
a full 2D — ► 3D crossover is expected in the spin flop. 
Notice that the MR is always negative, Ap/p < 0. 

In the Intermediate large distances regime, 
eo/(2Zt±) » K Q r ^> 1, the integral in |(5J) can be 
evaluated by direct expansion in powers of KorZtj_/eQ. 
Only the layers n = 0, ±1 contribute in this case, 
P(r,tx) = E n =o,±i iV'nl 2 , leading to 



P{r,t ± ) = g e - 2 ^l + 4( Ko r) 2 
for e /(2Zt±) > > 1 



(2*±) 5 



(11) 



From (JTTJ we obtain (r — > i?-r) across the spin flop 



&H>H 



f 



2/3 



(12) 



This formula corresponds to a crossover from a pure 2D 
case to an "intermediate dimension" (three transverse 
channels), and it is justified when {(JH>H f / &H<H f ) — ]- <C 
1. We have also performed a full numerical evaluation of 
P(r,t±) and of the MR with ip n (r) from Eq. ©, since 
the above considerations are based on asymptotic behav- 
ior. The exact numerical form was used in both Figures 
[Hand [31 below. We have found quite clearly that indeed 
in the temperature range where experimental data are 
available, T > 10K, the intermediate asymptotic formula 
(fl"2)) is the relevant one, since the MR is still relatively 
small there, e.g. \Ap/p\ m 0.35, T = 10K. However 
as the temperature is lowered, T < 10K, the MR in- 
creases and the system enters a crossover region between 
the asymptotic formulas Eq. (p~2|) and Eq. (fT0|) . which 
requires the use of the exact wave-functions. The cal- 
culated MR is plotted in Fig. [2] and we observe that the 
fitted value Ztx/eo = 0.057 agrees quite well with the es- 
timate Zt±/eo ~ 0.03 from band structure calculations^ 
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FIG. 3: (Color online.) The MR p(H)/p{0) (solid, blue line) 
for in-plane field H\\b, with Zt±/eo = 0.047. The squares are 
data from Ref. [U Inset: Evolution of the angle (|13|) with 
field. 



V. MAGNETORESISTANCE FOR IN-PLANE 
FIELD H || b 

In this case, due to the alignment of the DM-induced 
moments with the magnetic field, the spins rotate in 
the b — c plane (in opposite directions on the two sub- 
lattices) ji2ii2ii2ii5iiSi22i and align completely along c at a 
field H C 2- Thus, once we introduce the angle 9(H) that 
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the spins make with the ^-direction, our previous calcu- 
lations leading to Eq. (fT2")) remain valid with the replace- 
ment ij_ — > t± s'md. Observe that while for the undopcd 
LCO an intermediate in-plane spin flop is expected at 
H c i < H c2 r^± in LSCO the doped holes contribute to 
enhance the 6-axis spin susceptibility and then to confine 
the spins in the b — c plane, as it has been discussed re- 
cently in Ref. HH. We can then write the field dependence 
of sind between H = and H = H C 2 as 



smB{H) 



HD/ao 



AL t +4r?-(l-:cx™ P )tf 2 ' 



(13) 



In the above expression D is the DM anisotropy, rj = 
2JJ± is the inter-layer exchange, A out is the out-of-planc 
(or XY) anisotropy gap, <jq is the staggered order param- 
eter, x is the doping and Ximp is a dimcnsionlcss mea- 
sure of the holes-induced spin susceptibility. Here H is 
measured in units of g^usH, with g b s = 2.1 and fig is 
the Bohr magneton. The parameter values can be ex- 
tracted from the experiments , 15 ' 22 and for x = 0.01 we 
take D = 2.16 meV, r) = 1 (meV) 2 , A out = 3.2 meV, 
Ximp = 80, and er = 0.36, which gives H c2 « 19 
Using Zt±/ea = 0.047 as the only fitting parameter, we 
find a remarkable agreement with the experimental data 
of Ref. [13 at T = 20 K, as shown in Fig. El 

We also note that in Oxygen-doped compounds the 
contribution of the localized holes to the longitudinal 
susceptibility is much smaller, with Ximp ~ Ij^ 2 - As a 
consequence, one expects to observe at a field H c \ ~ 10 
T an intermediate flop which reflects in a kink in both 
the sm 6(H) curve and in the MR curve, as it has been 
measured indeed in the earlier transport measurements 
in La2Cu04+ y >i£ 



VI. DISCUSSION AND CONCLUSIONS 

In conclusion, we have shown that in the strongly lo- 
calized VRH regime for T < 50K, the in-plane MR is 



sensitive to the inter-layer hopping t± because an exter- 
nal magnetic field effectively changes the dimensionality 
of the problem, making the hopping quasi-2D. The MR 
reflects the physics of the spin-flop and is always negative; 
its value as well as temperature and field dependence are 
in excellent quantitative agreement with recent experi- 
ments in LSCO at x = 0.01 . 11 ' 12 Orbital effects, typically 
causing positive MR^ are negligible at this small dop- 
ing because the hole's localization length 1/k ~ 2 — 3 is 
much smaller than the magnetic length at any reasonable 
fields. 

Finally we comment that unlike LSCO where the role 
of magnetic anisotropies is very well established, in insu- 
lating YBa 2 Cu30 6+:r (YBCO) such anisotropies are ex- 
pected to be very weak, which is related to the absence of 
strong localization in this material^ Spin-related effects 
in the in-plane MR are absent for field in the c direction, 
while the in-plane field MR remains small and appears 
to be due to the dynamics of holes that are very weakly 
influenced by disorder ^ The complete understanding of 
these phenomena in YBCO remains an open issue al- 
though it is clear that magnctotransport is not domi- 
nated by the local spin physics. 
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